Transverse NMR Relaxation as a Probe of Mesoscopic Structure 
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Transverse NMR relaxation in a macroscopic sample is shown to be extremely sensitive to the 
structure of mesoscopic magnetic susceptibility variations. Such a sensitivity is proposed as a novel 
kind of contrast in the NMR measurements. For suspensions of arbitrary shaped paramagnetic 
objects, the transverse relaxation is found in the case of a small dephasing effect of an individual 
object. Strong relaxation rate dependence on the objects' shape agrees with experiments on whole 
blood. Demonstrated structure sensitivity is a generic effect that arises in NMR relaxation in porous 
media, biological systems, as well as in kinetics of diffusion limited reactions. 



NMR as a structure probe is utilized in the fields as 
diverse as chemistry, materials science, geology and 
biomedicine. Pristine specimen found in nature, such as 
rocks or biological tissues, possess complex structure at 
a mesocopic scale. This structure is of primary interest 
in numerous applications. For example, rock porosity in 
geology is important to assess the oil basin quality. In 
biological tissues the mesoscopic scale is set by the size of 
cells and blood vessels whose properties carry significant 
diagnostic and physiological information. 

It is the NMR monitored diffusion that is commonly 
accepted as a probe of mesoscopic structure in both in- 
organic [1] and living [2, 3] specimen. In the present 
Letter we propose a magnetic susceptibility contrast as a 
structure probe. Susceptibility inhomogeneities are often 
connected to the geometric structure, such as pore walls 
in porous media. In biological tissues they are brought 
by paramagnetic cells, such as deoxygenated red blood 
cells (RBCs) and iron-enriched cells in the brain gray 
matter. In some cases the susceptibility contrast can be 
artificially manipulated. 

Below we consider NMR signal from a suspension of 
arbitrarily shaped weakly paramagnetic objects. We 
demonstrate a significant individual object shape depen- 
dence of the transverse relaxation rate. We discuss this 
result in the biomedical context. Applications of the 
biomedical NMR imaging (MRI) are limited by a spatial 
resolution ~ 1 mm, which is larger than the cell size by 
2-3 orders of magnitude. Direct resolution enhancement 
is unfeasible since today MRI hardware hits physiologi- 
cal limits. Our results suggest that further progress can 
be made by a deeper analysis of the NMR signal since it 
contains significant information about the paramagnetic 
tissue structure at the scale of several /im. 

We compare our results with experiments on whole 
blood [4-6], with the objects being paramagnetic RBCs. 
Previous theoretical efforts in this context were focused 
on the effect of paramagnetic inclusions of specific geome- 
tries (spheres [7, 8] or cylinders simulating blood vessels 
[7, 9, 10]). The effect of object shape was not studied the- 
oretically although experiments [11, 12] and the Monte 
Carlo simulation [12] indicate a strong shape dependence 



of the transverse relaxation. 

We model the medium by a suspension of N 3> 1 identi- 
cal mesoscopic paramagnetic objects which are randomly 
placed and oriented. The NMR signal is acquired from 
nuclear spins that freely diffuse in the solvent and in the 
objects. A macroscopic volume V of the suspension is 
characterized by the volume fraction £ = Nvq/V of ob- 
jects (vq = single object's volume). The case of different 
object species is easily accounted for when (< 1, since 
they contribute additively to the relaxation rate [7, 9]. 

Transverse relaxation occurs due to two different mech- 
anisms: (i) microscopic spin-spin interactions at the 
molecular level, (ii) diffusion of spins in the magnetic 
field induced by mesoscopic objects. Fast processes (i) 
average out to produce a monoexponential relaxation. 
Processes (ii) are described in terms of the transverse 
magnetization density ^>(r), which evolves due to molec- 
ular diffusion and spin precession with the local Larmor 
frequency varying in space. It obeys the Bloch-Torrey 
equation [13] 

^= (W-^-^- iW (r)^, (1) 

where D is the diffusion coefficient of the molecules that 
carry spins, and Ti is the relaxation time due to the mi- 
croscopic interactions. The relaxation rate I/T2 is insen- 
sitive to magnetic field inhomogeneities at the mesoscopic 
scale. Rather, it characterizes local chemical composi- 
tion. We assume that D and T 2 are the same inside the 
objects and in the solvent. The constant term ujl pro- 
vides the Larmor precession in the homogeneous main 
field, and aj(r) = Y^l=i w o( r — r n) is the deviation from 
u>l due to the local magnetic fields induced by randomly 
located paramagnetic objects (as described later). 

The signal S(t) from a macroscopic sample is the sum 
of all spin magnetic moments regardless of their initial 
positions and their Brownian trajectories after the excita- 
tion [9]. In terms of the Green's function ip(r, To,t) of (1), 
defined by the initial condition ip(ro, r, t = Q) = 6(r— Tq), 

S(t) = l[d 3 vd 3 r V(r,r ,i) = e -***-*/r a s{t) . (2 ) 
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Microscopic processes decouple due to Eq. (1): -0 
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effects. The corresponding signal attenuation factor 
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describing these effects is the main object of our focus. 

Consider the mesoscopic part Af(r ,i) = 
fd 3 r 0(r o , r, t) of the spin packet magnetization. 
The d 3 r integration in (3) effectively averages M(r ,i) 
over randomly positioned objects. For £ <C 1, M is a 
product of factors contributed by individual objects [9]. 
In this case s(i) is expressible in terms of a single object 
dephasing effect f(t) [7, 9]: 



s = e-<f, f(t) 

y v o \ j / 1 

(4) 

Here rj is the mesoscopic part of the spin packet magne- 
tization density in the presence of a single object, 

^ = (^V 2 -^ (r))7?, 7 ? (r ,r,t=0)=<5(r-r ), (5) 

and average ( ) in (4) is taken over the object's orienta- 
tions. 

In the main field Bqz each paramagnetic object induces 
a local Larmor frequency shift uio(r) that is determined 
by the object's susceptibility profile x( r )- Below we use 
uniformly magnetized objects to compare with experi- 
ments: x(r) = x-v(r), x ^ 1> where v(r) is a shape 
function: v — 1 inside and v = outside the object. A 
convolution in r, lu in the Fourier space is 



Wo(k) = Stu ■ Y(k) ■ u(k) , Slu = Attx^l 



(6) 



where K(k) = 1/3 — k 2 z /k 2 is the longitudinal projection 
of an elementary magnetic dipole field, and the object's 
form factor w(k) is the Fourier transform of v(r). 

Transverse relaxation is qualitatively different in the 
limits of strong and weak dephasing. Introduce effective 
object radius p as that of a sphere of a volume vq. Water 
molecules pass by the object during the diffusion time 
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Typical phase acquired by the spins is 5u> -to- In the pre- 
sent work we focus on a weak dephasing case Sco-tD *C 1 
(diffusion narrowing regime). This regime covers a vari- 
ety of experiments, in particular spin dephasing in dia- 
magnetic and paramagnetic samples in the field Bq< IT. 

We find the Green's function n of Eq. (5) perturba- 
tively in the small parameter 5uj ■ tp, and use Eq. (4) 
to obtain f(t). This approach is analogous to the Born 
series for the quantum mechanical scattering amplitude. 

The zeroth order in 8lo ■ to describes free diffusion. In 
this case the total magnetization of each spin packet is 



conserved, J d 3 r r?(ro, r; t) = 1 in (4), and s(t) — 1. The 
first order correction to / vanishes since it is proportional 
to the angular average of the dipole field. The expression 
for / is dominated by the second order in 8u> - to (Fig. 1): 
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Here r is the dimensionless time r = i/to. The inner 
integral in Eq. (8), which is taken over the directions of 
q = kp, depends exclusively on the object shape. The 
object size enters Eq. (8) only through the diffusion time 
tj3, Eq. (7). The function g depends on the particular 
sequence of the radiofrcquency (rf) pulses applied to ma- 
nipulate the spins and is discussed later. 

As a conservative estimate, the formal series for fit") 
converges when a < 1. Indeed, each successive term in 
the perturbative expansion of f(t) is multiplied by the di- 
mensionless factor ia — i47r/(27r) 3 (6ui ■ tn) (vo/p 3 )- An- 
gular integrations of the products Y(ki)Y(k.2)... improve 
convergence by bringing additional factors < 1 that are 
object shape specific and hard to account for in general. 

Odd orders of the expansion in a are imaginary. They 
renormalize the homogeneous component of the suspen- 
sion's magnetic susceptibility. Since the first order van- 
ishes the correction to u>l is proportional to a 3 . The 
signal attenuation is determined by the even orders in a. 
The correction to (8) is of the order of a 4 and is negative. 

Consider the free induction decay (FID), an evolution 
after a single rf tt/2 pulse which creates the maximal 
transverse spin magnetization. The function g in (8), de- 
noted as <7fiDj is proportional to a time convolution of the 
three free diffusion propagators (q, t) = 6 , (r)e -9 r , 
9{t) being a unit step function (Fig. 1, left): 
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To reduce sensitivity to large scale field inhomo- 
geneities, samples are often irradiated by a number of re- 
focusing rf 7r pulses. Each such pulse quickly rotates the 
spins by 7r around an axis which is transverse to z. This 




FIG. 1: Second order processes for fit). Left: FID relaxation. 
Circles, wavy lines and crosses stand for —iSui, ^(k), and v(k) 
respectively. Solid lines represent free propagators ry'°'(k, r) 
in time intervals between interactions. External momenta are 
set to zero due to Eq. (3). Right: CPMG relaxation. Each 
section represents a free propagator r/^ in the interval At 
between successive refocusing pulses. Complex conjugation 
on every other interval At is indicated with the filled circle. 
Equation (12) is obtained as a sum of all such configurations. 
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is equivalent to a complex conjugation of rj developed up 
until this time moment. The resulting distribution rj* is 
the initial condition for the further evolution. 

In the spin echo (SE) technique [14] a single tt pulse is 
applied at the time tg/2 and the signal is measured at 
t = tg. The corresponding g function reads 
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In the CPMG (Carr-Purcell-Meiboom-Gill) protocol [15] 
refocusing n pulses are generated in a long train and the 
steady state signal is studied as a function of the inter- 
pulse interval At (Fig. 1, right): 



ffCPMG = q t - 2 tanh 
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Equations (10)-(12) yield that at r < 1, / oc r 2 for 
the FID and / oc r 3 for the SE and the CPMG sequences. 
Asymptotic expansion of (8) in t- 1 / 2 at t > 1 gives 
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with t4 F id = Vtt, ^se = (2^2 - l)\fn, and Aqpmg = 
(2\/2- 1) C(3/2) /Vtt ~ 2.695 for the considered pulse se- 
quences. The dimensionless NMR relaxivity r-2 is shown 
in Fig. 2, left, for the case of the homogeneously mag- 
netized spherical particles. Shape dependence is illus- 
trated in Fig. 2, right, for the case of disk-shaped objects. 
The height-to-radius ratio c defines the disk shape, with 
c = 0.5 being close to the intact RBC. 

Below we analyze our results, Eqs. (8) and (13). 

(i) Relaxation (4) crucially depends on the shape of 
the object. It is the form factor 6(k) that governs the 
convergence of the integral for large q = kp in (8). The 
integral converges at k ~ 1/p, allowing one to probe the 
object's structure. (A quantum mechanical analogy is 
scattering amplitude dependence on the form factor of 
the external potential.) A point- like magnetization v oc 
6(r) causes a divergence in Eqs. (8) and (13). In the 
present case this "nonrenormalizability" (non-universal 
cutoff dependence) effectively increases sensitivity in the 
NMR measurements. 

(ii) Shape sensitivity is a consequence of a singular 
interaction Y ~ r~ v between nuclear spins and objects. 
Consider the case when the singularity in Y is cut off at 
a scale r < a. Then Y(k) — » as ka > 1. If a > p, the 
integral in Eq. (8) is insensitive to the form factor since it 
converges at k<l/a<l/p, destroying shape sensitivity. 
Physically, such a cutoff introduces a spherical "cloud" 
of a radius a around each object. This cloud smears in- 
formation about the object's structure. The power nec- 
essary for shape dependence is v > 2 + (d — 2)/N for the 
N-th order in d dimensions. Thus both magnetic dipole 
(y = 3) and contact interaction Y = 6(r) in d — 3 yield 
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FIG. 2: Mesoscopic relaxivity T2 = /(t)/t for a 2 = 15/2-7T. 
Left: Objects are spheres, t = tB for the SE, t — At for the 
CPMG. Right: Shape effect: disks vs. spheres. CPMG relax- 
ivity r2(At/to = 1), objects are disks with height-to-radius 
ratio c, and spheres of the same volume. 



shape sensitivity already in the second order, as shown 
above. 

(iii) Shape sensitivity is present for any held Bq- 
Above we demonstrated shape sensitivity in the domain 
where the perturbative approach is reliable (a<l). We 
now prove it for any a. Integrals such as (8) whose con- 
vergence is form factor dependent appear in each order 
of the perturbation series for f(t). Although angular in- 
tegrations impede explicit summation of this series, they 
do not cause nonanalyticity at a = 0, and thus radius of 
convergence in aotBo is finite. Therefore the series can 
be analytically continued to the large field domain a > 1 
where the perturbation theory formally breaks up. The 
final result for f(t) would still be form factor dependent, 
which proves shape sensitivity for any field. 

(iv) Shape sensitivity is a generic effect.- Consider a 
diffusion limited chemical reaction on impurities with a 
shape u(r). The FID signal analytically continued by 
—iuj(r) — > u(r) gives the impurity shape dependent reac- 
tion rate. 

Below we compare the results (8) and (13) with exper- 
iments. As a test we use the reported relaxation rate in 
dilute (C = 0.02) suspensions of polystyrene microspheres 
in paramagnetically doped water [16] (Fig. 3, left). 

Further experiments were performed on the deoxy- 
genated blood with a high RBC volume fraction £ = 
0.40 - 0.60 [4-6]. To apply Eqs. (8) and (13) one needs 
to take into account a slower diffusion inside the cells 
and to extend our approach for large £. The former will 
be considered elsewhere. For now, we obtain upper and 
lower estimates for the relaxation rate by using in Eq. (8) 
the values Di n , D ont of the diffusion coefficient in erythro- 
cytes and plasma respectively. 

The £ ~ 1 case poses a challenging task equivalent to 
finding the statistical sum of a dense gas of objects. In- 
stead we replace £ by £(1 — £) in (4), which is well sup- 
ported experimentally [17]. Such a replacement is justi- 
fied by the virial expansion. Eq. (4) treats exactly the 
first cumulant of the statistical sum. The second cumu- 
lant provides an C(£ 2 ) negative correction. This together 
with a vanishing mesoscopic contribution as £— >1 justi- 
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FIG. 3: Theory (lines) vs. experiments (symbols). Left: Re- 
laxation rate ife = —(\ns)/t for the FID (boxes) and SE (cir- 
cles) as a function of the particle diameter. Filled and hollow 
symbols correspond to measured and Monte Carlo simulated 
relaxation rates [16] . Right: CPMG relaxation rate for the hu- 
man blood samples [6] for B = 1.41, 1.18, 0.94, 0.71 T (from 
top to bottom). Experimental errors are 10-20% [6]. Follow- 
ing our discussion after Eq. (8), theory agrees with experiment 
for small a and overestimates it for a ~ 1. 



lies a quadratic polynomial interpolation, ((1 — (). The 
latter is correct as C^0 and 1, and describes a crossover 
between the dilute and the extremely dense cases. 

The relaxation rate in deoxygenated blood measured 
in [4] quadratically depends on the magnetic field, 
— (\ns)/t = ki-Bq, in agreement with Eq. (8). The pro- 
portionality coefficient n\ was found to be 7.2 s _1 T~ 2 for 
the CPMG pulse sequence with At = 4 ms. In [4], the 
field range B a = 0.05 - 1.5 T yields a = 0.033 - 0.99. 
We calculated ( — 0.55 from the parameters given in 
[4], utilized the magnetic susceptibility of the deoxy- 
genated RBCs x = 2.7 x 10~ 7 [5], and simulated the 
intact erythrocytes by disks of the known volume of 
87 /zm 3 with the height-to-radius ratio of c = 0.5. Using 
D ant = 2.20 /im 2 /ms and An = 0.76 jUm 2 /ms [18], our 
theory gives 4.7 < Kith < 5.6 s _1 T~ 2 . 

To assess this result we note that neither the suscep- 
tibility of RBCs nor their actual shape was reported in 
[4]. Chemicals used to treat the samples are likely to 
change osmotic pressure in plasma, which would deform 
the RBCs thus changing all relevant parameters. 

Quadratic dependence of the SE blood relaxation rate 
on x, which follows from Eq. (8), was confirmed by vary- 
ing the RBC oxygen saturation y in the field B — 1.5 T 
[5]: — (\ns)/t — k 2 (0.95 — J/) 2 , with the measured co- 
efficient K2 = 55 s _1 for C = 0.3 and re 2 = 59 s _1 
for £ — 0.4. Our approach results in the corresponding 
ranges 26 < K 2tn < 56 s _1 and 30 < K 2t h < 64 s -1 . 

The CPMG relaxation rate R 2 = - (In s)/i in the 
whole blood was measured [6] as a function of the in- 
terecho interval (Fig 3, right) for 0.71< J B <1.41 T. We 
simulated blood as described above using D — D out for 
plasma. The use of the value D m instead of D out yields 



about the same rate i? 2 for the short times and approxi- 
mately a twofold increase of i? 2 for the large times. 

This brief survey shows that, although crude, our 
model captures essential features of the NMR relaxation. 
Experiments at higher fields [11, 12] confirm the shape 
dependence for a > 1. Their results can be well de- 
scribed by adjusting to and a [8] or by fitting to a sim- 
ple chemical exchange model [19]. However, fitting has a 
predictive power when the signal universally depends on 
a handful of phcnomcnological parameters. Shape sen- 
sitivity makes such a fitting meaningless in the case of 
varying tissue structure. Because of the same reason, in 
experiments analogous to [4-6, 11, 12] it is essential to 
control volume fraction, shape and susceptibility of para- 
magnetic objects, and effective diffusion coefficient in the 
sample. 

To conclude, we showed that transverse relaxation 
from a suspension of paramagnetic objects is extremely 
sensitive to the shape of the individual object. This sen- 
sitivity to geometric structure is a generic effect that can 
be employed as a novel type of contrast in NMR mea- 
surements thus effectively increasing spatial resolution. 
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